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Abstract 

We present a preliminary measurement of the inclusive charmless semileptonic branching fraction 
of the B meson, based on 140 fb~^ of data collected by the Belle detector at the KEKB e+e" 
asymmetric collider. Events are tagged by fully reconstructing one of the B mesons, produced 
in pairs from T(4S'). The signal for h u semileptonic decay is distinguished from the h ^ c 
semileptonic background using the hadronic and leptonic invariant mass distributions Mx and q^. 
We find the partial branching fraction for the kinematical region given by Mx < 1-7 GeV/c^ and 
^2 > 8 GeVVc2, Ai3(5 ^ XJv) = [0.99 ± 0.15(stat) ± 0.18(syst) ± 0.04(6 ^ u) ± 0.07(6 ^ 
c)] X 10~^. Using a theoretical prediction for the extrapolation to the full range of Mx and 
variables, we obtain B{B XJu) = [3.37 ± 0.50(stat) ± 0.60(syst) ± 0.14(6 u) ± 0.24(6 
c) ib0.50(/„ error)] x 10^^. From this measurement, the magnitude of the the Cabibbo-Kobayashi- 
Maskawa matrix element Vub is |Kfe| = [5.54 ± 0.42(stat) ± 0.50(syst) ± 0.12(6 ^ u) ± 0.19(6 ^ 
c) ± 0.42(/„ error) ± 0.27(^ ^ error)] x 10"=^. 

PACS numbers: 12.15.Hh,11.30.Er,13.25.Hw 
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INTRODUCTION 



Measurement of the Cabibbo-Kobayashi-Maskawa matrix element | Vub \ is crucial to test 
the Standard Model prediction oi CP violation. In this paper, we report a preliminary 
measurement of the inclusive charmless semileptonic branching fraction of the B meson 
B{B Xuiu), which provides one of the best ways to determine \Vub\- In this measurement, 
one of the B mesons, referred to as the tag side meson, -Btag; is fully reconstructed in several 
hadronic decay modes to tag the production, flavor and charge as well as the momentum of 
the B meson. The semileptonic decay of the other B meson, referred to as the signal side 
meson, Sgig, is detected by the presence of a high momentum electron or muon. This method 
allows one to reconstruct the invariant mass of the hadronic system in the semileptonic decay, 
Mx- The invariant mass squared of the leptonic system, g^, can also be determined by 
inferring the missing neutrino momentum. Both kinematical quantities are used to separate 
the B Xu^v signal decays from the abundant B — > Xclf background decays, with a good 
signal-to-noisc ratio. A similar type of analysis was performed by BaBar, where only a cut 
on Mx was applied [1]. According to a recent theoretical calculation [2], a simultaneous 
cut on the two variables is beneficial because it reduces theoretical uncertainties. Such a 
simultaneous cut on the two variables was used in a similar analysis by Belle where, instead 
of the full reconstruction tagging, a simulated annealing technique was applied to separate 
the two B meson decays [3]. 

This measurement has an advantage over the measurements based on the lepton momen- 
tum alone [4] in that it covers a relatively large phase space of the signal spectrum, which 
reduces the error in extrapolating the measured partial branching fraction to obtain the 
total branching fraction. While the method has these unique features, it requires a large 
sample of BB events because the full reconstruction efficiency is rather low, typically of the 
order of 0.3%. 

The data used in the present analysis were collected with the Belle detector at the asym- 
metric energy KEKB [5] storage ring. The Belle detector [6] is a largc-solid-angle magnetic 
spectrometer that consists of a three- layer silicon vertex detector (SVD), a 50- layer central 
drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like 
arrangement of time-of-fiight scintillation counters (TOF) , and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a super-conducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instru- 
mented to detect mesons and to identify muons (KLM). 

The result presented in this paper is based on a 140 fb~^ data sample recorded at the 
T(45') resonance, which contains 152 x IQ^BB pairs. An additional 15 fb~^ data sample 
taken at a center-of-mass energy 60 McV below the T(4S') resonance is used to subtract the 
background from the e+e^ qq process (g = m, rf, s, c). 

Monte Carlo simulated events were used to determine the efficiency as well as signal 
and background distributions in the control variables. The simulation program is based on 
GEANT [7] and fully describes the detector geometry and response. To model the B X^iv 
decays, we employ a combination of exclusive channels, where X„ is either a vr or a p [8] 
or an excited Xu state [9], and an inclusive model for non-resonant final states [10]. The 
B Xciiy transitions are simulated according to the QQ decay generator [11]. 

This paper is organized as follows. We first discuss the reconstruction of the tag side 
B meson, requirements on the signal side and the signal yield extraction. The observed 
number oi B ^ X^ii^ decays is then converted into a partial branching fraction for the 
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kinematical region given by Mx < 1-7 GeV/c^ and > 8 GeV^/c^, for which the theoretical 
uncertainties in modehng h ^ u transitions are small. The theoretical result [2] quoted above 
together with the 6-quark shape function parameters determined from the Belle 6 — > S7 
measurement [12] is used to extract the branching fraction for charmless semileptonic decays, 
which is, finally, used to obtain the magnitude of the Vuh matrix element. 

RECONSTRUCTION OF THE TAGGING SIDE 

In the present analysis, i?tag candidates are reconstructed in the decay modes, — >■ 
D^*^-TT+/p+/at/Dl*^+ and B+ D^*^^7t+ / p+ /at /Dl*^+ . Inclusion of charge conjugate 
decays is implied throughout this paper. 

Primary charged tracks are reconstructed with hit information from the CDC. They 
are required to satisfy track quality cuts based on their impact parameters relative to the 
measured profile of the interaction point of the two beams. 

Charged kaons are identified combining specific ionization (dE/dx) measurements in the 
CDC, Cerenkov light yields in the ACC and time- of- flight measurements in the TOP. The 
kaon identification efficiency is approximately 88% and the average pion fake rate is about 
8%. 

Candidate 7r° mesons are reconstructed using 77 pairs with an invariant mass between 
117.8 and 150.2 MeV/c^. Each photon is required to have a minimum energy deposit of 
> 50MeV {E^ > 30McV for neutral pions from D* decays). 

Kg mesons are reconstructed using pairs of charged tracks that have an invariant mass 
within ±30MeV/c^ of the known Kg mass and a well reconstructed vertex that is displaced 
from the interaction point. Candidate p"^ and p° mesons are reconstructed in the decay 
modes tt'^tt^ and 7r"^7r~, by requiring their invariant mass to be within ±225MeV/c^ of the 
nominal p mass. Then, af candidates are selected by combining a p° candidate and a pion, 
if their invariant mass lies between 0.7 and 1.6GeV/c^ and if the three tracks form a good 
vertex. 

r>° candidates are reconstructed in the — > K^n', ir+7r~7r°, K'^n'^ 
KgTr'^Ti^, /C^7r+7r"7r'^ and K^K~ decay modes, while D~ candidates are reconstructed in the 
D~ K^n^n^ , K^n^Tr'TT^, K^n" , fC°7r~7r°, KgiT'ir'n'^ and K^K~ti~ decays. candi- 
dates are reconstructed in the decay modes — > KgK^ and K^K'n^ . These candidates 
are required to have an invariant mass mo within ±4 — Su of the nominal D mass, where the 
mass resolution a depends on the decay mode. D* mesons are reconstructed by combining 
the D candidate and a low momentum pion, D*~ D^n^ /D^tt^ and D*^ Z)°7r°. They 
are required to have a mass difference Am = mj^^ — within ±5 MeV/c^ (±4 — Gcr) of 
the nominal value. Por the decays with a photon, D*^ — > D^^ and D*"*" — > -D^7, we require 
that the mass difference Am = rriQ^ — rriQ be within ±20 MeV/c^ (=1=2(7) of the nominal 
value. 

The selection of Btag candidates is based on the beam-constrained mass, Mbc = 

^Je^LJc^ -Pb/c^ and the energy difference, AE = E% - E^^^. Here E4„„, = ~ 
5.290 GeV is the beam energy in the center of mass system, and p*^ and E*^ are the cms 
momentum and energy of the reconstructed B meson. (Throughout this paper the variables 
calculated in the center of mass system will be denoted with an asterisk.) Events satisfying 
Mbc > 5.22GeV/c^ and \AE\ < 0.3 GeV are subject to further analysis. 

The combinatorial background from jet-like e'^e~ — > qq processes is suppressed by event 
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FIG. 1: Beam-constrained mass (Mbc) distribution for the B~^{a) and B^{h) candidates with 
the — 0.2GeV< AE < 0.05 GeV requirement. The scaled distribution of off-resonance data is 
subtracted. The shaded areas indicate the results of the fit in the Mbc signal region. 

topology cuts based on the normalized second Fox- Wolfram moment {R2) [13], R2 < 0.5, 
and for some modes also by a cut on | cos6'thrust| < 0.8, where 6'thrust is the angle between 
the thrust axis of the B candidate and that of the rest of the event. 

The signal region for the tagging B is defined with the cuts Mbc > 5.26GeV/c^ and 
— 0.2GeV< AE < 0.05 GeV. The cuts are optimized to maximize the statistics of the 
signal, while minimizing the migration of tracks between the tag and signal sides. 

If an event has multiple Bt^g candidates, we choose the candidate having the smallest 
based on the deviations from the nominal values of AE, mo-, and Am if applicable. 

Figure 1 shows the distribution of Mbc for the B^ and B^ candidates in the AE signal 
region. The qq background contribution is subtracted using the scaled off-resonance data. 
The numbers of tagged events are estimated to be (1.58 ± 0.18) x 10^ for B^ and (2.47 ± 
0.22) X 10^ for 5+, by fitting the distribution with empirical signal [14] and background 
functions [15]. The reconstruction efficiencies are 0.21% and 0.33%, while the purities are 
47% and 50% for the B^ and B^ samples, as determined using scaled off-resonance data 
subtraction. 



RECONSTRUCTION OF THE SIGNAL SIDE 

For events tagged by fully reconstructed i?tag candidates, we search for electrons or muons 
from semileptonic decays of the signal side B meson. Electron identification is based on a 
combination of the dE/dx value as measured in the CDC, the response of the ACC, the 
shower shape in the ECL and the ratio of the energy deposit in the ECL to the momentum 
measured by the tracking system. Muon identification by the KLM relies on the number 
of hits in resistive plate counters interspersed in the iron yoke. The lepton identification 
efficiencies are about 90% for both electrons and muons in the momentum region above 
IGeV/c. The hadron fake rate is measured using — > vr+vr" and (j)^ — *• K^K~ decays, 
and found to be less than 0.2% for electrons and 1.5% for muons in the same momentum 
region. 
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FIG. 2: Beam-constrained mass (Mbc) distribution for prompt semileptonic decays (a) and events 
satisfying all B Xu£iy signal requirements, including Mx < 1.7 GeV/c^ and > 8 GeV^/c^ 
(b). The shaded areas indicate the results of the fit in the Mbc signal region. 

We select electrons having p* > 0.6GeV/c and pt > 0.6GeV/c, and muons with p* > 
0.8GeV/c and pt > 0.7GeV/c, with pt being the transverse momentum component with 
respect to the positron beam axis. We also require that they are detected in the laboratory 
polar angular region of 26° < 9 < 140°. Backgrounds from J/ip decays, photon conversions in 
the material of the detector and vr*^ Dalitz decays are minimized by imposing veto conditions; 
we calculate invariant masses for each lepton candidate when combined with opposite charge 
leptons and with an additional photon in the case of electrons {rriee-y), and reject the 

lepton if rriu hes within ±5(7 of the nominal J/ip mass, rrige within ±100MeV/c^ or rriee-y 
within ±3cr of the nominal vr*^ mass. 

For the prompt semileptonic decay signal, we require a lepton with momentum p* exceed- 
ing 1 GeV/c. We also require the lepton charge to be consistent with a prompt semileptonic 
decay, when the Stag candidate is charged. No requirement is imposed on the lepton charge 
when the Stag candidate is neutral. 

The B Xuiv signal events are selected by imposing several requirements to suppress 
the B Xciv background as well as the background from poorly reconstructed events. Fake 
charged tracks, arising mainly from duplication in the tracking of low momentum curling 
tracks, are removed on the basis of the angle between the two track candidates, and the 
difference in their momenta. We require that the event contains only one lepton, has zero 
net charge Qi = 0) and has a missing mass consistent with zero (—1 GeV^/c^ < m^j^^ < 
0.5GeV^/c^). In order to further suppress the B — > Xf.ii' background, we require that the 
number of kaons, either K"^ or Kg, is zero {Nk = 0) on the signal side. To reject events 
containing a meson, we require that the angle between the missing momentum and the 
direction of the candidate cluster be greater than 37°. 

Finally, we select the signal events with requirements on the invariant mass of the hadronic 
system Mx and invariant mass squared of the leptonic system g^. The variable Mx is 
calculated from the measured momenta of all charged tracks and energy deposits of all 
neutral clusters in the ECL that are not used in the Stag reconstruction or not identified 
as leptons. Charged tracks are assigned the mass of the pion, kaon or proton, based on the 
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FIG. 3: Mx distribution for the selected events with > 8 GeV^/c^, data and fitted contributions 
of the B Xclv and B X^ii^ transitions. 

information from the particle identification system. According to Monte Carlo simulations 
the resolution in Mx for the selected events is found to be about 125 MeV/ and 130 MeV/ 
for B —>■ Xud-v and B —>■ Xciiy processes, respectively. The four-momentum transfer q is 
calculated as g = Pt(4S) — PBtag ~ Px, where Pt(4S), PStag Px are the four-momentum 
vectors of the T(45'), -Stag and the reconstructed hadronic system, respectively. 

Our signal region is defined as Mx < 1-7 GeV/c^ and > 8 GeV^/c^. With these 
requirements, we suppress the kinematic region where the theoretical interpretation of the 
result would be difficult [2], and also reduce the level of the B —>■ Xciiy background. The 
distribution in the beam constrained mass for the selected events is shown in Fig. 2(b). 

SIGNAL YIELD EXTRACTION 

The signal yield is extracted by a fit to the distribution of hadronic invariant mass Mx 
for the selected events with > 8 GeV^/c^. The distribution of events in hadronic mass is 
determined by dividing the data into several Mx bins. For each bin, the yield is extracted 
from a fit with empirical signal [14] and background functions [15] to the corresponding Mbc 
distribution, an example of which is shown in Fig. 2(b). Figure 3 shows the resulting Mx 
distribution. 

To determine the raw number of events corresponding to the B X^^i^ process, we 
first fit the Mx distribution with two contributions, the signal and the background from 
the B — > XJm process. The shapes of both contributions are determined from Monte 
Carlo simulation. The result of the fit is indicated in Fig. 3. Using the same relative 
normalization of the processes, we plot in Figs. 4 and 5 the distributions of events over the 
Mx (finer granularity as compared to Fig. 3) and variables. In the regions of low Mx 
(< 1.7GeV/c^) and high q^ a clear contribution from B X„£z/ can be observed. 

The raw yield of 5 — > X„£z/ decay in the signal region, A*"^^ = 174 ± 26, is determined 
by subtracting the fitted B — ^ Xciv contribution (Fig. 3). The error is statistical only. 

In order to extract the partial branching fraction AB for B — >■ in the signal region 

Mx < 1.7 GeV/c^, > 8 GeV^/c^, Monte Carlo simulation is used to convert the observed 



8 




FIG. 4: Mx distributions for the selected events with > 8 GeV^/c^: a) data and fitted contri- 
butions of the B Xclv and B Xuii' transitions, b) measured distribution for the B Xuiv 
transition, obtained by subtracting the fitted B XcH-v contribution from the data, compared to 
the simulated B Xuiv contribution with the normalization given by the fit. The final cut on 
Mx is indicated by a dotted line. 




FIG. 5: The q'^ distribution for the selected events with Mx < 1.7 GeV/c^: a) data and contribu- 
tions of the B — > Xcii' and B Xu^v transitions, normalized with the fit to the Mx distribution, 
b) measured distribution for the B X^i-v transition, obtained by subtracting the B X^iv 
contribution from the data, compared to the simulated B — > Xu(-v contribution. The final cut on 
is indicated by a dotted line. 



number of events Nl^^ to the true number of signal events produced in this region, Nf,^u, 
and to estimate the efficiency for these events to be observed anywhere. Nb^u is estimated by 
Nb^u = A^fe!f„ X F, with F = 1 + N2/N1 - N3/N1. Here A^i is the number of simulated events 
observed in the signal region and N2 {N^) is the number of events generated inside (outside) 
the signal region and observed outside (inside) the signal region. We find F = 0.984 ±0.014, 
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and thus N^^u — 171 ± 26. The efficiency e^^" for selecting B — > Xulv events after the 
lepton momentum cut is predicted to be 27.4%. 

The relative partial branching fraction AB{B — > Xuii') / B{B X£i>) is obtained by 

ABjB^XJu) 1 etc 4 

B{B X£u) 7V,i eLt" e^-" ef-" " ^ ^ 

Here A^si — (5.07 ± 0.04) x 10^ is the number of events having at least one lepton with 
p* > l.OGeV/c, determined from a fit to the corresponding Mbc distribution (Fig. 2(a)), 
and corrected for the expected fraction of background events from non-semileptonic decays 
(14.1%), as estimated by MC simulation. The factor efl^c/^frec accounts for a possible dif- 
ference in the ^tag reconstruction efficiency in the presence of a semileptonic or 5 — > X^ii^ 
decay; ef/e\~*'^ is the ratio of fractions of semileptonic decay leptons with p* > 1 GeV/c, in 
the whole phase space for the B XJ,v events, and within the region Mx < 1.7 GeV/c^ 
and > 8 GeV^/c^ for signal events. The product of efficiency ratios is found to be 
4eo/4Z^ X /e^" = 0.75 ± 0.048. 



SYSTEMATIC ERRORS 

The major sources of systematic error in the relative branching fraction are the uncer- 
tainty in the background subtraction to extract the yield N^^u, the uncertainties in the 
calculation of the efficiency ratios, and the uncertainty due to the treatment of 5 — >■ X^ii^ 
decays in the Monte Carlo simulation that is used to estimate e^^^ and F. 

The validity of the B — > Xciiy background simulation is tested with the B — > Xciu 
enhanced control sample, where all selection cuts are applied with one exception: we require 
at least one kaon in the event. We have checked that for the Mx distribution of this control 
sample there is good agreement between data and the simulation. The relative errors arising 
from uncertainties in the background subtraction to extract A^^^u are estimated to be 7% 
due to the uncertainty in the modeling of the B — > Xf.ii' background and 15% due to the 
limited statistics for the simulated background events. 

Statistical uncertainty in the determination of efficiency ratios contributes 6% to the 
systematic error. The uncertainty in the ratio F/ Cg^f " is estimated by varying the parameters 
of the simulation model for B — > Xuii^ decays. We assign a 4% error after varying the 
parameters of the inclusive model within their errors. 

The uncertainties due to systematic errors in tracking efficiency, particle identification 
efficiency and cluster finding efficiency are estimated by varying them by their respective 
errors, and observing the effect on the value of partial branching fraction. These sources 
give correlated errors on simulated b ^ u and b c events so we add or subtract the error 
of each source linearly for the two samples, according to the relative sign of their effect. The 
contributions for different sources are combined in quadrature to a systematic error of 6.5%. 
The total systematic error excluding b ^ u and ft — > c model dependences amounts to 18%. 
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SUMMARY 



In summary, the relative partial branching fraction for B — > Xuii^ decays in the kinematic 
region Mx < 1.7 GeV/c^, q'^ > 8 GeV^/c^ is 



[0.92±0.14(stat)±0.17(syst)±0.04(6^ 'u)±0.06(6^ c)]xlO-^ , (2) 



where the first error is statistical, the second experimental systematic, and the third and 
fourth arise from uncertainties in modeling the B Xuii^ and B Xcii' transitions. By 
using the measured semileptonic branching fraction B{B — > Xiu) — 0.1073 ± 0.0028 [16], 
we obtain 

AB{B^XJiy) = [0.99±0.15(stat)±0.18(syst) ±0.04(6 ^m) ±0.07(6 ^c)]x 10-^ . (3) 

The branching fraction B{B — > Xuii') is calculated from the above value through the 
expression 

B{B ^ XJu) = AB{B ^ XJu)/U ■ (4) 

The extrapolation coefficient /„ is estimated to be 0.303± 0.035 using the De Fazio and Neu- 
bert (FN) prescription [10] with the 6-quark shape function parameters = 4.62 GeV/c^ 
and /x^ = 0.40 GeV^/c^ and their one sigma error ellipse, which are determined from a 
recent Belle B — > X^'y measurement [12, 17]. It is rescaled to 0.294 ± 0.035 by multiplying 
by a factor /„o(BLL)//„o(FN) = 0.324/0.334, where /„o(BLL) and /„o(FN) denote the /„ 
values calculated for the values of = 4.71 GeV/c^ and /x^ = 0.2 GeV^/c^ by the Bauer, 
Ligeti and Luke (BLL) [2] and Dc Fazio and Neubert prescriptions, respectively. In both 
as and l/m(, expansions the BLL prescription contains corrections of higher order than the 
FN prescription. The uncertainty of is modified by including the contributions from the 
subleading shape function and the weak annihilation. The former is estimated to be 4% [18] 
and the latter is estimated to be 8% [2]. As a result we use /„ = 0.294 ± 0.044 in Eq.(4). 
This yields 

B(B XJv) = [3.37 ± 0.50(stat) ± 0.60(syst) ± 0.14(6 ^ u) (5) 
±0.24(6 ^ c) ± 0.50(/„ error)] x 10"^ , 

where the additional error comes from the uncertainty in the calculation of /„. 
Combining this result with the average B lifetime Tb = (1.587 ± 0.011) ps [16], 
the CKM matrix element \Vub\ is obtained by using the formula 

|V;.|= 0.00424 (B(B^XJ.)1M^Y' 

' V 0-002 Tb J ^ ^ 

which is an updated version of the expression given in [19], and incorporates the latest 
measurements of the heavy-quark parameters [20]. Finally, we obtain 

\Vub\ = [5.54 ± 0.42(stat) ± 0.50(syst) ± 0.12(6 ^ u) ± 0.19(6 ^ c) (7) 
±0.42(/„ error) ± 0.27(i3 ^ |Kb| error)] x 10"^ . 

The first four errors are statistical, systematic, b ^ c and b ^ u model dependence. The 
latter two are due to uncertainties of /„ and of the relation between B and |Kb| in Eq. (6), 
respectively. 
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The present work demonstrates the effectiveness of Mx and measurements of 5 — > 
Xud-v decays using full reconstruction tagging. The result is consistent with previous mea- 
surements of B{B — > Xu^l') [1, 3, 4, 21]. In the future, further accumulation of data as well 
as MC data will allow us to reduce the statistical error, and better understanding of the 
signal and the background components will help to improve the experimental systematics 
as well as constrain the theoretical uncertainties. 
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